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INTRODUCTION 


The  physical  principles  of  the  deterministic  properties  of  surface 
waves  have  been  studied  for  over  a century.  The  work  by  Lambf is  one  of 
the  classic  studies  in  this  fields  The  consideration  of  the  statistical 
properties  of  the  ocean  surface  was  first  started  about  20  years  ago. 

One  of  the  first  publications  on  this  subject  was  by  Pierson,  Neumann 
and  James2.  This  manual  is  informative  and  readable.  ^Two  excellent  ^ 
surveys  of  ocean  statistics  have  been  given  by  KinsmanJ  and  Cartwright  . 
One  of  the  original  and  much  quoted  work  on  surface  statistics  was  done 
under  project  SWOP  in  1960^.  A conference  in  ocean  wave  spectra0  contains 
the  work  nf  many  sMpntiBt.fi  in  the  field. 

The  purpose  o?Hhis_study  6s— fe«Preview^fche  work  on  ocean  surface 
statistics^*  In  this" review  some  new  approaches  are  considered.  A 
brief  discussion  of  some  similarities  between  ocean  surface  statistics 
and  underwater  sound  noise  fields  is  included. It  is  hoped  that  this 
review  will  be  helpful  as  an  introduction  to  the\field. 
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Single  Frequency  Plane  Waves 


Let  us  first  consider  a single  frequency  plane  wave  moving  in  the 
direction  . Let  the  amplitude  of  the  wave  be  Q.* > and  the  phase 
€n  • Let  'KM,#,*')  be  the  displacement  of  the  water  height  from 
its  equilibrium  position,  where  % and  ft  represent  a point  on  the  surface 
and  £ is  the  time.  Then 


In  this  equation,  we  have  assumed  that  the  surface  wave  is  for  the  deep 
water  case  and  for  this  case  the  wave  number  , where 

ft  is  the  gravity  acceleration  constant.  Equation  #r  defines  the  sur- 
face height  for  all  points  on  the  surface  and  for  all  time.  It  is  a 
three  dimensional  wave,  two  of  the  dimensions  are  the  x and  y coordinates 
and  the  third  dimension  is  the  displacement  of  the  surface  in  the  V 
direction  given  by>| 6).  Lei  us  now  assume  that  there  are 
single  frequency  plane  waves  of  the  same  frequency  each  travelling  in 
a different  direction.  Then  the  displacement  is 


In  equation  2,  we  have  assumed  that  the  principle  of  superposition  holds 
i.e.  that  the  total  displacement  is  equal  to  the  linear  sum  of  the 
individual  displacements.  In  general  this  is  not  true.  However,  for 
small  displacements,  the  approximation  of  assuming  linear  superposition 
is  good.  In  addition,  it  makes  much  of  the  theory  considerably  simpler 
and  many  experimental  results  may  be  predicted  from  the  linear  super- 
position model.  A given  plane  wave  travelling  in  the  &•»  direction,  has 
an  amplitude  and  phased  . Up  to  this  point  the  description  has 
been  deterministic.  We  now  introduce  the  random  model  by  considering 
an  ensemble  of  surfaces  similar  to  the  deterministic  surface.  Each 
surface  consists  of  /V  single  frequency  plane  waves  *he  fr**  plane 
wave  is  travelling  in  the  direction  with  an  amplitude  <h%  . For  the 
/V**  plane  wave  in  a given  surface,  the  phased,  is  selected  from  a 
population  with  phase  between  O and  Hit*,  all  phases  being  equally 
likely.  In  this  model,  each  surface  is  identical  except  for  the  phase 
(rn  • IF  we  obtain  the  probability  distribution  of  the  phase  over  an 
ensemble  of  surfaces,  then  the  probability  density  is  shown  in  Fig.  1. 
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Fig.  1. 

Because  of  the  random  phase  condition,  each  surface  displacement  config- 
uration is  different  than  the  others.  It  should  be  noted  that  the  model 
could  be  constructed  by  having  random  amplitudes  Cbn  in  addition  to 
random  phase,  or  random  amplitudes  alone. 

Let  us  now  obtains  ome  of  the  statistical  properties  of  the  ensem- 
ble. However,  before  we  do  this,  let  us  review  some  definitions  that 
will  be  used.  A stationary  random  procer.s  is  one  whose  statistical 
properties  do  not  change  with  time.  An  ergodic  process  is  a stationary 
process  such  that  with  probability  one,  the  time  average  of  any  single 
function  of  the  ensemble  is  equal  to  the  ensemble  average.  It  should 
be  noted  that  a process  to  be  ergodic  is  that  each  sample  of  the  ensem- 
ble go  through  all  possible  points  of  the  set.  The  model  as  described 
above  is  stationary  but  not  ergodic.  For  example,  a given  surface  is  a 
superposition  of  sinusoidal  plane  waves  of  a given  frequency,  so  that 
it  is  a sinusoidally  oscillating  surface  in  time  of  the  same  frequency 
as  its  components.  The  configuration  that  the  individual  surface  can 
take  is  quite  broad,  but  the  frequency  of  each  individual  spatial  point 
is  still  sinusoidal.  However,  the  amplitude  of  the  oscillations  at  the 
point  is  fixed  and  will  thus  not  necessarily  reach  all  possible  values. 
There  are  several  ways  to  circumvent  this  problem.  However,  the  method 
of  single  frequencies  is  very  helpful  intuitively  and  since  the  process 
will  become  ergodic  when  continuum  of  frequencies  is  introduced,  we  will 
use  it  here. 

He  will  now  average  over  the  ensemble  of  systems.  In  order  to  hold 
the  amount  of  symbols  to  a minimum,  we  will  consider  a two  dimensional 
wave  surface.  The  extension  to  a three  dimensional  surface  is  obvious, 
once  the  two  dimensional  surface  is  understood.  For  a two  dimensional 
surface, 

tUXlt)^'ZanOy<,(,i^X-co-6- €*) 

0*!  d 
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We  now  assume  that  the  process  is  stationary  in  the  vide  sense  in  both 
space  and  time.  That  is  the  second  order  properties,  such  as  correla- 
tion, depend  only  on  the  difference  in  x and  the  difference  in  t and 
not  on  where  we  begin.  We  let  yj  r „ %.  - £■  _ A 


- ' Hhnllo  Cut  - Co-0  ( — 

The  symbol^*  ) denotes  the  expected  value  over  the  ensemble.  We 
now  use  the  trigonometric  identity  that  relates  the  product  of  cosines 
to  the  sum  and  differences  of  the  arguments  of  cosines.  Then 

hlx-Ojt-ioy  = -Krxa-~a,cosCy  x*  f-Cx-o-vr- 

We  note  that  the  first  term  on  the  right  hand  side  is  always  zero  when 
averaged  over  € . That  is,  the  phase  difference  remains  when/»n  jiff 

and  when  » n,  < c«CO.*6^  - O 

When#**»^/)  the  second  term  is  also  zero,  but  when  the  phases  in  the 
second  term  cancel,  we  have 

#13#  ^ V 

The  mean  square  surface  height  & is  defined  as 


We  now  define  the  spatial  temporal  correlation  of  the  surface  height  as 

Ho.  v.  -ri  * <" -it*.  vj9-yj  *-r)> 


Then  using  similar  arguments  as  in  the  two  dimensional  case 

Equation  4,  related  the  spatial- temporal  correlation  of  the  ocean  sur- 
face, for  single  frequency  waves,  to  the  value  of  which  is  propor- 
tional to  the  power  of  a given  component  wave  in  the  d*  direction. 

We  can  now  extend  this  approach  to  include  other  frequencies.  Then 
d,  is  a function  of  direction  0 and  frequency  tO  . Finally,  we 
can  let  this  double  sum  increase  and  let  the  increments  between  values 
decrease  so  that  the  double  sum  increase  and  let  the  increments  between 
values  decrease,  so  that  the  double  sum  approaches  the  double  integral. 
An  approach  of  this  type  has  been  given  by  Pierson  and  is  reported  by 
Kinsman^.  Kinsman's  equation  in  our  notation  is 
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called  the  directional  wave  spectrum.  Ax O^j  • ) ^ & 

is  the  relative  amount  of  power  of  the  component  surface  waves  at  C ^ d) 
in  the  range  A three  dimensional  plot  of  /\*-(  U>,  6)  is  given 

by  KinsmazH  in  Pig.  8.  3-1.  This  equation  5 represents  the  relation 

between  the  spatial- temporal  surface  correlation  and  the  directional 
wave  spectrum. 

In  this  section,  we  have  gone  to  great  lengths  to  treat  the  single 
frequency  case  and  to  obtain  the  relationship  given  by  equation  4.  We 
have  indicated  that  the  same  type  of  reasoning  could  be  used  to  obtain 
equation  5- 

Another  form  of  equation  5 that  is  useful  is  to  express  the  right 
hand  side  in  terms  of  wave  numbers  rather  than  and  & . 

Let  A 


Then 


we  obtain  the  following  two  equations 


If  we  use  these  transformations 
The  first  relation  is 


The  second  equation  is 
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Properties  at  a Single  Point 

is  very  difficult  to  obtain  experimentally.  The  method 
of  attack  as  explained  in  the  literature,  is  to  consider  various  sim- 
plifications of  this  quantity.  That  is,  various  parameters  are  held 
fixed  during  a given  measurement. 

The  simplest  quantity  to  measure  is  the  time  variation  of  the  sur- 
face at  a given  point.  We  use  an  instrument  to  measure  the  height  ver- 
sus time  at  a given  point  of  the  ocean  surface.  A capacitance  wave  pole 
is  an  example  of  this  type  of  instrument.  We'  then  obtain  a time  history 
as  shown  in  Fig.  2. 


From  this  time  history,  we  obtain  the  autocorrelation  as 

o/>0  ^ <0)  Co.  o.  O A C °j  °j  4- >0> 

The  power  spectrum  at  tnis  point  is  the  Fourier  transform  of  the  auto- 
correlation function.  Thus 

$(f) e Ji e*rC-ja irfr) 

A great  deal  of  the  work  in  ocean  statistics  and  the  literature  are 
concerned  with  the  quantity  aw-  Returning  to  equation  5>  we  have 

I jZo  /3u,  AVcw,  &)  eo+CvT*) 

eco, 0,1* ) 

Note  that  is  a double  sided  spectrum  and  is  single  sided. 

Now  consider  the  right  hand  side  of  the  last  equation,  i.e. 

AH**)  ~ AXC*>,  ©) 


'0C903 


Utl 


d)  is  the  directional  wave  spectrum.  In  the  equation,  we  inte 
grate  over  all  d . Thus  all  knowledge  of  the  direction  of  the  waves 
is  lost,  if  we  consider  only  the  quantity  • For  the  Neumann- 

Pierson  spectrum 

a*/.,  ^ _ c -zr  *e  4.  m #7 


Then 


= 0,  otherwise 
g s 980  cm/sec^  2 
U = wind  speed  in  cm/ sec. 

C » 3.  05x10*  cm2/sec.  5 
fjj  * is  a cutoff  frequency 
^ « o for  a fully  aroused  sea 
we  obtain  a curve  as  shown  in  Fig.  3. 


For  a higher  value  of  wind  speed  U,  the  peak  of  this  curve  shifts  to  a 
lower  value  of  w and  becomes  narrower.  If  we  define  © as 


where  fp  is  the  peak  frequency  and  Tc  and 
are  the  frequencies  at  the  half  power  points,  then  &■=./,*/  for  all 
values  of  U.  Thus  the  Neumann- Pierson  spectrum  represents  a wide 
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h 


band  process.  The  integration  over  in  equation  8,  gives  the  total 
power  at  a point.  Thus 


* f cm  exp(~  * 


owei- 

Properties  for  a Given  Time 


co 


Let  us  now  return  to  equation  #6  and  consider  the  spatial  corre- 
lation of  the  surface  height  at  an  instant.  Then 

<r'eo,  v$  o)  fcflL,  o v) 


We  can  now  express  } in  terms  of  its  even  and  odd 

components.  Thus  ® 


k c-4s  h~) 


Thus 


k - /yWc-', k o) co^(Xoh-H,v)iojy 

and 

fc  (A*,  A^~h-ff6^cu>y'T^  <^+ V>wrft; 


O',  >0 


Now  the  instantaneous  height  of  the  surface  may  be  obtained  by  photo- 
grammetric  techniques  as  shown  by  Pierson^.  From  the  instantaneous 
heights,  the  correlation  6 ( ^ itfj^nay  be  computed.  If  we  then  use 
equation  9,  we  may  numerically  integrate  the  right  hand  sides  of  this 
equation  and  thus  obtain  A*  and  A^  (•&.*,  ^3)  • Finally 

CJkn,  ie  obtained  from  its  even  and  odd  components.  It 

should  be  noted  that  ClV»*So)  may  be  obtained  by  stereo  methods. 
However.  requires  photographs  at  different  times.  This 

value  is  difficult  to  obtain  in  practice.  The  difficulty  is  overcome,  if 
it  is  assumed  that  the  directional  wave  spectrum  occurs  only  in  two 
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adjacent  quadrants  about  the  wind  axis.  That  is,  the  surface  wave 
directions  are  always  within  - ^ of  the  wind  axis.  Methods  of  this 
type  have  been  used  in  project  SWOP  (Stereo  Wave  Observation 
Project)®. 

For  analytic  purposes,  the  problem  is  simplified  if  it  is  assumed 
that  there  is  180°  symmetry  in  the  directional  wave  spectrum,  i.  e. 

. bL4.,,Jk») 

c'eio^o)--  fUJLJA}  no 

From  equation  10.  it  is  seen  that  the  two  dimensional  Fourier  trans- 
form of  A*'('&hjA9') is  Thus  A'£^4j^i8  the  inverse  two 

, dimensional  Fourier  transform  of  o)  That  is, 

'k’CJLcl£i ) - +*l9v)]  m 

i 

We  now  make  the  additional  assumption  that  the  surface  is  isotropic. 

An  assumption  of  this  type  was  made  by  W.  Marsh  and  also  by  Nuttall 
and  Cron®.  By  an  isotropic  surface,  we  mean  that  the  directional  wave 
spectrum  is  the  same  in  all  directions. 

Let  i r . nr.  7Z 


/r-  -s  J o*-+Y 


°(  ~ 4-4! 


Then  CC 

Let  A / J A j Ai  A)  and  O ■=  /r-  c*-*  & 

W's  & ogiir 

*)  ~ y/rx  / / ec++<t>  4-4*K«Ad 

oo  a ir 

Ai  640)  * Je  - <t>)J 

If  we  use  the  relation  for  the  Bessel  function  of  the  first  kind,  i.  e. 

rS  ■f  e*p(j*Ce~a'>  * 6 

Then  (.<<)  - /n- 

This  is  a Fourier- Bessel  transform.  The  inverse  transform  is  of  the 

same  form.  That  is  M 

•©  ^ 

- *7r /<H  *K  3IC4*-)  ^3 
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In  the  preceding  discussion  we  have  considered  the  rectangular 
form  of  the  directional  wave  spectrum,  i.  e.  Wethen 

found  the  relation  between  the  spatial  correlation  of  the  surface  and 
the  directional  wave  spectrum,  for  the  isotropic  case.  The  two 
quantities  are  Fourier-Bessel  transforms  of  one  another.  The  direc- 
tional wave  spectrum  is  usually  expressed  in  the  polar  form.  For 
example  the  Neumann-Pierson  spectrum  in  polar  form  is  given  in 
equations  #7  and  #8.  We  will  now  derive  the  relation  between  the  spa- 
tial correlation  and  the  directional  wave  spectrum,  expressed  in  the 
polar  form,  for  the  isotropic  case.  Let  T''—  o , in  equation  #5.  Then 

C'CCq K ei)* i P* /Vk <0 coo c c-*e  + ^ * ) 

-n  a o 9 


For  the  isotropic  case,  At'Cu^  A 

Liet  /r~  = i|  u«-  + v •*-  jT's.t • The 

fcr-)  X -i-  © /dw  | 

Then  using  the  definition  of  the  Bessel  1 
obtain  *o 

Ai/O  LA-  T* 


Equation  14  has  been  used  in  Marsh^  et  al.  The  inverse  Fourier- 
Bessel  transform  of  (/r'^  is 

Cron  and  Streit^  have  numerically  evaluated  equation  14,  for  the  case 
of  the  Neumann- Pierson  spectrum.  This  was  done  for  various  values 
of  wind  speed.  Nuttall  and  Cron®  have  numerically  evaluated  equation 
15  for  various  assumed  values  of  spatial  correlation. 


Let  us  now  consider  two  points  of  the  surface.  We  let  time  vary 
and  obtain  Q(0i  Vj  for  these  two  points.  The  cross  spectral 
density,  denoted  b^  ) in  this  article,  is 

9 < T)  rrfi  7-)  #16 

Note  that  U and  V are  held  fixed  in  the  integration  process.  Barber 
has  shown,  as  referenced  by  Cartwright*  that  a linear  array  of  ele- 
ments and  their  associated  cross  spectral  densities  can  be  used  to 
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obtain  the  directional  wave  spectrum.  For  example,  for  two  points  in 
space,  we  could  obtain  the  height  versus  time  of  each  point,  as  given 
in  Fig.  4. 


Fig.  4. 


^et  h C O,  o,  4) 

Then  r 

g'-pc "sKr)*  li~.  X j6 

T “9  "6 

Thus  Vj  T*)  can  be  obtained  experimentally.  The  Fourier  transform 

of  would  yield  Vj  f ) . An  alternate  method  is  to  FFT 

and  (fc)  so  as  to  obtain  N,ff)  and  ) respectively.  Then 

Vt  Wi*£)  where  represents  the  complex  conjugate. 

For  a large  number  oi  elements,  this  procedure  is  repeated  for  all 
pairs  of  elements. 

We  will  now  review  Barber's  method.  Let  be  the  power 

spectrum  at  a point.  The  autocorrelation  function  is 

- <hce.  o,fc )•)£<>, 

Thus  f>C©,o,T')  < > A^(f)  where  ^ * 

represents  the  Fourier  Transform.  We  again  ^ represent 
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the  directional  wave  spectrum  for  a plane  wave  of  this  power 
approaching  the  two  elements  as  shown  in  Fig.  5,  we  have 


The  right  hand  side  represents  the  cross  spectral  density  of  a plane 
wave.  We  now  consider  the  ocean  surface  to  consist  of  plane  waves 
from  all  directions.  Then  the  total  cross  spectrum  is 

- fAxCf,(»  expCij 4c/c^<9)o/e 

Since  Al(fj  is  periodic  in  & with  period  7.TT  , we  may 
express  it  as  a Fourier  series.  Then 


In  this  expression,  we  have  assumed  symmetry  in  6 . Physically, 
this  would  correspond  to  taking  the  wind  direction  as  the  axis  and 
measuring  © from  this  axis.  We  also  have  the  line  of  the  linear 
array  along  the  wind  axis.  For  the  general  case  we  need  two  linear 
arrays,  one  along  the  axis  and  one  perpendicular  to  the  axis.  Substi- 
tuting equation  17  into  equation  16  and  using  the  definition  of  the 
Bessel  function  which  is  xi r 

~ tf  j d S 

XTT  o 

We  then  obtain  , 

gC'tjvJ')  •=  jt  Ao  ©0  T0&O  + 37r  2 4 LAi 

We  can  measure  K'f')-  ^e  can  compute 

If  we  terminate  the  infinite  series  at  the  term,  we  have  N 

unknowns  in  A*U^')  where  i = o,  J, N - 1 . If  we  now  use 

N receivers,  we  will  have  N equations  and  N unknowns.  We  can  then 
solve  these  equations,  to  find  the  N values  of  A i Ifc)  . From  these 
N values,  we  can  compute  the  approximate  value  of  A\f,  O ) from 
equation  17.  Thus  cross  spectral  measurements  may  be  used  to  obtain 
the  directional  wave  spectrum. 


VML, 
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A variation  of  the  above  techniques  was  used  in  the  underwater 
sound  case,  by  Cron,  Hassell  and  Keltonic*®.  In  that  reference, 
various  directional  wave  fields  were  considered  and  corresponding  to 
each  wave  field,  the  spatial  correlation  of  a vertical  array  of  elements 
was  derived  analytically.  Cross  correlation  measurements  were  then 
made  on  a vertical  array  of  elements  and  the  experimental  values  were 
compared  with  the  values  obtained  from  the  theoretical  model  of 
various  assumed  directional  wave  spectra.  The  directional  wave 
spectra  that  gave  the  best  agreement  was  chosen  as  the  actual  direc- 
tional wave  spectrum.  Another  variation  of  this  technique  used  in 
underwater  sound  is  that  by  Von  Winkle  et  al**. 

We  will  now  consider  the  probability  distribution  of  the  ocean 
surface  height.  However,  before  we  do,  let  us  summarise  the  last 
section.  We  started  by  indicating  how  the  relation  between  the  spatial 
temporal  correlation  and  the  wave  directional  spectrum  could  be 
obtained.  The  spatial  temporal  correlation  was  expressed  as 

We  then  took  the  special  case  of  one  point  on  the  surface.  For 
this  case,  we  obtained  the  expression  This  expression  was 

related  to  the  power  spectrum.  The  directional  properties  of  the 
waves  were  lost.  We  then  specialized  the  general  spatial  temporal 
correlation  to  the  case  of  a fixed  instant  of  time  and  thus  obtained 
QCVj  Vj  oj  . This  expression  was  related  to  the  wave  directional 
spectrum.  We  then  considered  two  fixed  points  on  the  surface  and 
obtained  (?(  °j  K,  ”*0.  F rom  C (w,  *0  ),  the  cross  spectral  density  was 

obtained. 

All  of  these  properties  are  second  order  properties.  For 
example,  tpis  a second  order  property.  We  have  not  covered 

higher  order  moments  such  as^VUi^eJ)  and^i1* «,  * 
which  are  third  order  moments.  If  the  process  is  Gaussian,  only  the 
first  and  second  order  moments  are  needed.  Fortunately,  for  the  sur- 
face height  distribution,  the  Gaussian  process  may  be  assumed  for 
many  cases. 


Probability  Distributions 


Let  us  now  consider  the  probability  distribution  of  the  surface 
height  at  a point.  We  will  return  to  the  principle  of  linear  superposition. 
The  height  of  the  surface  at  a point  is  the  sum  of  many  sine  waves 
travelling  in  different  directions  and  with  different  frequencies.  For  a 
broad  set  of  conditions,  the  sum  of  a large  number  of  random  vari- 
ables results  in  a random  variable  with  a Gaussian  distribution,  as 
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may  be  obtained  from  the  central  limit  theorem.  The  probability  die 
tribution  of  the  height  is  v 

where  h is  the  height  and  G*-*  <V>  is  the  mean  square  height.  This 
distribution  has  a bell  shape  and  is  shown  in  Fig.  6. 


Fig.  6. 

It  should  be  noted  that  the  principle  of  superposition  holds  for  small 
amplitude,  low  frequency  waves.  High  frequency  waves  result  in  non 
linear  interaction  terms  and  thus  one  would  not  expect  a Gaussian 
process.  This  is  borne  out  by  experiment.  ' 

Mac  Kay,  as  referenced  by  Kinsman3,  made  an  exhaustive  study 
of  the  wave  height  distribution,  by  the  means  of  bottom  pressure 
recorders.  He  used  the  chi-square  test,  Kolmogorov-Smirnov  test 
etc. , to  show  that  the  Gaussian  hypothesis  was  accepted.  However, 
it  is  known  that  the  high  frequency  surface  waves  are  greatly  attenu- 
ated. Thus  the  high  frequency  surface  waves  were  not  included  in 
MacKay's  experiments.  Kinsman3  obtained  the  probability  distribution 
of  the  surface  height  by  means  of  a capacitance  pole.  He  thus  included 
the  high  frequency  components.  A comparison  between  the  experimen- 
tal values  and  the  Gaussian  distribution  showed  that  the  two  were 
different,  but  only  slightly  different.  There  is  a small  amount  of 
skewness,  i.  e.  The  third  moment  was  not  zero  in  the  experimental 
values,  whereas  the  Gaussian  distribution  is  symmetric.  Since  the 
two  distributions  are  close,  for  some  practical  purposes,  we  can 
assume  that  the  surface  height  will  also  be  Gaussian. 

If  the  surface  is  stationary  in  the  wide  sense  in  both  space  and 
time  and  if  the  process  is  Gaussian,  then  the  directional  wave  spectrum 
or  the  spatial-temporal  surface  correlation  completely  specifies  the 
statistical  properties  of  the  surface. 


V.  ___  SUMMARY 

0ufJiiiSS 

We  have  outlined^the  method  for  obtaining  the  equation  showing 
the  relation  between  the  spatial-temporal  correlation  and  the  direc- 
tional wave  spectrum.  Methods  of  measuring  the  spatial- temporal 
correlation  of  the  surface-ws^^Vnentloned  along  with  the  equations,  a* 
Special  cases  o/jjj&e  spatial-tempo  rancor  relation  wePetconsidered. 
The  most  important  one^Irom  the  present  measurement  viewpoint (i.  e 
the  simplest  to  measure)  is  the  temporal  correlation ^pt  a point  on  the 
surface.  The  Fourier  Transform  of  this  correlatioi 


leads  to  the 

power  spectrum  of  the  Neumann- Pierson  type.  Finally,  a method  of 
obtaining  cross  spectra  between  pairs  of  a linear  array  of  elements  '-s 
waF  discussed  and  similarities  between  this  method  and  the  under- 


water sound  problem  of  obtaining  the  directionality  of  ambient  noise 


—was  mentioned. 
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